Diquark and Pion Condensation in Random Matrix Models for two-color QCD 
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We introduce a random matrix model with the symmetries of QCD with two colors at nonzero 
isospin and baryon chemical potentials and temperature. We analyze its phase diagram and find 
phases with condensation of pion and diquark states in addition to the phases with spontaneously 
broken chiral symmetries. In the limit of small chemical potentials and quark masses, we reproduce 
the mean field results obtained from chiral Lagrangians. As in the case of QCD with three colors, 
the presence of two chemical potentials breaks the flavor symmetry and leads to phases that are 
characterized by different behaviors of the chiral condensates for each flavor. In particular, the 
phase diagram we obtain is similar to QCD with three colors and three flavors of quarks of equal 
masses at zero baryon chemical potential and nonzero isospin and strange chemical potentials. A 
tricritical point of the superfluid transitions found in lattice calculations and from an analysis in 
terms of chiral Lagrangians does not appear in the random matrix model. Remarkably, at fixed 
isospin chemical potential, for the regions outside of the superfluid phases, the phase diagram in the 
temperature - baryon chemical potential plane for two colors and three colors are qualitatively the 
same. 



I. INTRODUCTION 

Recently, the phase diagram of QCD at nonzero chem- 
ical potential for the baryon number has received a great 
deal of attention. It has been found that a plethora of 
phase transitions are possible as a function of the baryon 
chemical potential, ^b, the isospin chemical potential, 
/i/, and the temperature T. Among others, a supercon- 
ducting phase 0,13, a pion condensateJJQ and a critical 
endpoint have been predicted P, H, Q (see for a re- 
view). A systematic analysis of the phase diagram is pos- 
sible in the fiB = plane and, perturbatively, for asymp- 
totically large values of Ql- At finite, not asymptoti- 
cally large baryon chemical potential, where the interac- 
tions between the quasi-particles that carry color charge 
are strong, the theoretical basis of these predictions is 
rather weak. So far, it has not been possible to con- 
firm the existence of the proposed phases in this region, 
neither from experimental results nor from lattice cal- 
culations. The reason for the failure of standard lattice 
simulations is that in the presence of a chemical poten- 
tial for baryons the fermion determinant in the partition 
function becomes complex so that a probabilistic inter- 
pretation of the partition function is not possible. 

Thus, it is attractive to try to gain insight into the 
properties of such condensed phases by studying QCD- 
like theories which avoid some of these complications. In 
this paper we study QCD with two colors with baryons 
that are colorless diquark states. As in the case of three 
colors, chiral symmetry is spontaneously broken and the 
theory is confining at low temperatures and chemical 
potential. However, at nonzero baryon number chemi- 
cal potential, two important simplifications occur in this 
theory. First, the condensing diquark states are Gold- 
stone bosons associated with the spontaneous breaking of 
chiral symmetry with a critical chemical potential equal 



to half the pion mass. This makes it possible to ana- 
lyze the condensed phase by means of chiral Lagrangians 
I3, 0,0,0,^]- Second, as a consequence of the pseu- 
doreality of the SU{2) gauge group, the fermion deter- 
minant remains real when a chemical potential is intro- 
duced. Therefore, for an even number of flavors, this 
theory can be simulated on the lattice without encoun- 
tering a sign problem. 

The phase diagram in two-color QCD has been stud- 
ied using a variety of methods. The appearance of a 
baryonic diquark condensate above a critical chemical 
potential oi ^bc — 'TfiTi/'i was first observed in lattice 
simulations at strong coupling 0, 0, 0|. More re- 
cently, the phase diagram was analyzed in the frame- 
work of chiral Lagrangians j^l, ilCV ilL il2i] . where it 
was found that the critical value of the chemical po- 
tential is indeed given by the mass of the lightest 
state, and that the transition is of second order at zero 
temperature. Lattice simulations confirm these results 
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In [23,123,133 it was observed that beyond a critical value 
of the baryon chemical potential and temperature, the 
transition to the diquark condensation transition ceases 
to be of second order and becomes a first order transi- 
tion. This implies the existence of a tricritical point. In 
chiral per turbation theory, such a tricritical point was 
found jTll [3^ at values for the chemical potential and 
the temperature of the order of the pion mass m,T^. 

In physical applications, such as in neutron stars or 
in heavy ion collisions, the chemical potentials for differ- 
ent quark flavors are not necessarily equal. While the 
chemical potential for baryon charge is the same for all 
quarks, an additional isospin chemical potential parame- 
terizes the difference between the chemical potentials for 
two light quark flavors. Its effect on the phase diagram 
has been studied for Nc = 2 in 13] using effective chi- 
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ral Lagrangians, and for iVp = 3 in 33] using a random 
matrix model, in |34ll35! . 36] by means of a Nambu-Jona- 
Lasinio model, and in '37] within the ladder approxima- 
tion. 

In this article, we employ a random matrix model to 
study QCD with two colors at finite chemical potentials 
for isospin and baryon density and finite temperature. 
Originally, random matrix models have been introduced 
in QCD in order to describe correlations of low-lying 
eigenvalues of the Dirac operator [s^ ■ These random 
matrix models are equivalent to the static part of a chiral 
Lagrangian which is uniquely determined by the symme- 
tries of the microscopic theory l4ll li^ . lisj ] . From this 
connection, it has been shown that they provide an ex- 
act analytical description of correlations in the low-lying 
Dirac spectrum. In addition to their use as exact ana- 
lytical models for the spectrum, random matrix models 
can also serve as schematic models for the phase transi- 
tions in QCD. The first of these models was introduced 
in [3 , fisf to describe the chiral restoration transition at 
finite temperature. At finite baryon chemical potential, 
they have been successfully used to explain the failure of 
the quenched approximation in lattice QCD ^46j , and to 
predict a tricritical point for the chiral restoration transi- 
tion at finite density and temperature . The static part 
of the effective chiral Lagrangian for QCD with the cou- 
pling of charged Goldstone bosons to a chemical potential 
has also been derived from a chiral random matrix model 
I47I ] . More recently, closely related models with a random 
gauge potential have been considered for the supercon- 
ducting phases |49l l50j ] and for diquark condensation for 
Nc = 2 |51|. Though these models have additional sym- 
metries beyond those considered in this article, they lead 
to the same spectral density as the corresponding chiral 
random matrix ensembles |53 ]. 

The random matrix theory approach rests on the idea 
that it is the matrix equivalent of a Landau- Ginzburg po- 
tential. Both theories are based on the symmetries of the 
microscopic partition function. Indeed, as we will also see 
in this paper, it is possible to derive a Landau-Ginzburg 
functional from the random matrix theory. The pres- 
ence of a diquark and a pion condensation phase raises 
several interesting questions, which can be addressed by 
studying a random matrix model. In particular, we hope 
to investigate the robustness of the tricritical point that 
has been found in lattice simulations |3Clj] and in chiral 
perturbation theory [T]| . In general, we wish to iden- 
tify which features of the phase diagram are generic and 
due to the underlying symmetries and which features are 
model dependent. 

This article is organized as follows. In section^ we 
give a short review of the symmetries and phase diagram 
of the partition function of QCD with two colors. Our 
random matrix model is introduced in section IIIII We 
show that it can be rewritten exactly in terms of the 
static part of a chiral Lagrangian. With an ansatz based 
on the structure of the expected condensates we obtain 
an effective potential for these condensates. Our main 



results are presented in section where we analyze the 
phase diagram in the chemical potential plane. 

Concluding remarks are made in section IVIl 

II. QCD WITH TWO COLORS 

The partition function for QCD with two colors is given 

by 



^QCD — ^J|det(Z? + m/ + /x/7o)^ , (2.1) 

where the Euclidean Dirac operator is D = "ifiD^^ = 
7/i(c^/i + 7m ^be Euclidean 7-matrices, and A^^ 

is the S'f7(2)-valued gauge potential. For each of the Nf 
quark flavors, the quark mass is given by m/ and the 
chemical potential by /z/. The brackets denote an aver- 
age over the Euclidean Yang-Mills action. In the case of 
only two flavors, the chemical potentials can be expressed 
in terms of isospin and baryon chemical potential, 

f^i = ^(^ii - Ai2)- (2.2) 

Due to the pseudoreahty of the SU{2) gauge group, 
the flavor symmetry group is enlarged to SU{2Nf) from 
the usual SU{Nf) x SU{Nf) in case of QCD with three 
or more colors. In the vacuum, this symmetry is spon- 
taneously broken to Sp{2Nf ). Thus, in the chiral limit, 
one expects Nf{2Nf — 1) — 1 exactly massless Goldstone 
bosons. Among these Goldstone bosons there are di- 
quark states which carry baryon charge, in addition to 
the mesonic pion states which carry a charge with re- 
spect to isospin number. The mass term breaks the 
chiral flavor symmetry explicitly down to Sp{2Nf). A 
baryon chemical potential term by itself is invariant un- 
der an U{1)b X SU{Nf ) x SU{Nf ) subgroup of the orig- 
inal SU{2Nf), and in the presence of a mass term this 
symmetry is further broken down to U{1)b x SU{Nf). 

The presence of an additional chemical potential cou- 
pling to light bosonic states greatly increases the number 
of possible phases. At low temperature and chemical po- 
tential we expect that chiral symmetry is spontaneously 
broken. The phases with broken chiral symmetry are 
characterized by the order parameters (uu) and (dd) . As 
in the case of A^c = 3 ,33] , in the presence of two different 
chemical potentials for the two flavors u and d, there is 
no reason to believe that (uu) — (dd) , and in general, one 
expects that the chiral condensates for both flavors are 
not equal. At /i/ ^ 0, when the two chemical potentials 
are different, it cannot be excluded that chiral symmetry 
restoration will take place via two separate phase transi- 
tions. 

In addition to the pion condensate p — ■^{{uj^d) — 
(d75u)), already known from the case Nc — 3, which 
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appears above a critical isospin chemical potential, for 
QCD with two colors we also have a diquark conden- 
sate ji{{d^u*) — {u^d*) + (d^u) — (u'^d)) above a criti- 
cal baryon chemical potential. From general arguments, 
one expects a phase transition when the chemical po- 
tential becomes equal to the ratio of the mass and the 
charge for the lightest charged excitation. Indeed, this 
has been shown explicitly by an analysis in terms of chiral 
Lagrangians at low chemical potential and temperature 
H In 111 III El- It was also found that the transi- 
tion to a baryonic diquark condensation phase or a pion 
condensation is of second order. 

From the symmetries of the partition function alone, it 
is possible to predict some of the properties of the phase 
diagram. As has been shown in |9, 10], the determinantal 
term in the partition function (|2.1|) can be rewritten in 
such a way that the enlarged flavor symmetry becomes 
manifest. Using a chiral representation for the Euclidean 
7-matrices, in which jojfj, = diag(i(T^, — icr^) with the 
spin matrices = {—i,ak), and the pseudoreality prop- 
erty of the generators of the gauge group, —T2TaT2 = t*, 
the fermionic part of the Lagrangian with the Dirac op- 
erator can be rewritten in terms of a new spinor basis. In 
the chiral limit to/ — *■ 0, but including a nonzero chemi- 
cal potential, it becomes. 







M/ 
icr,,D,, 



(2.3) 



where the spinors = <y2T2'4^*ji transform in the same 
representation of the color group as the spinors ipL- 
Clearly, at /iy = and to/ — 0, the flavor symme- 
try group is enlarged to SU{2Nf). For two flavors, the 
fermion determinant in the partition function factorizes 
into 

det(TO^ - {iDfj,a^ + Hi)a2T2{i{Df_,a^)'^ - fii)<J2T2) 

X det(TO^ - {iDf_,a^ + fi2)<y2T2{iiD^a^)'^ - H2)<J2T2)- 

(2.4) 

The partition function is therefore invariant under the 
interchanges Hf —>■ — /x/, / = 1,2, and /ii ^ 112- Trans- 
lated into the physical chemical potentials, the partition 
function is invariant under the individual sign changes 
fii — i- —fij and fiB ~^J■B, and under the interchange 
fJ-i ^ fJ-B- The latter symmetry, which amounts to 
II2 — > — /i2, means the interchange of quarks and conju- 
gate anti-quarks. Under this transformation, a diquark 
condensate is transformed into a pion condensate. Thus, 
overall one expects that the phase diagram is symmetric 
about the lines /ij = fi^ and the lines = and iib = 0- 
Therefore from the symmetry of the QCD partition func- 
tion we can restrict our study of the phase diagram to 
one half quadrant in the /ib — /x/-plane. 

III. RANDOM MATRIX MODEL 

We will introduce a random matrix model for QCD 
with two colors and study its phase diagram at finite 



chemical potentials for isospin and baryon number and 
finite temperature. The idea of the random matrix model 
is to replace the matrix elements of the Dirac operator by 
Gaussian distributed random variables, while respecting 
all the global symmetries of the QCD partition function. 
The dependence on temperature and chemical potentials 
will enter through external fields, which arc also subject 
to symmetry constraints. 

Besides the chiral and flavor symmetries, the Dirac op- 
erator can also have additional anti- unitary symmetries. 
For an SU{2) gauge group, the Dirac operator satisfies 
the commutation relation 



[iD, CT2K] = 0, {CT2Kf - 1, 



(3.1) 



where C = 7274 is the charge conjugation matrix, T2 is a 
generator of the gauge group, and K is the complex con- 
jugation operator. Due to this anti-unitary symmetry, it 
is possible to find a basis in which all matrix elements 
of D are real. Thus, the matrices in the corresponding 
random matrix ensemble will have real entries, which is 
denoted by the Dyson index (3=1. 

In the absence of a temperature and chemical poten- 
tial, the matrix elements of the Dirac operator are simply 
replaced by real random variables with a Gaussian dis- 
tribution. The chemical potential breaks the flavor sym- 
metry in the same way as in the QCD partition function. 
Just as the dependence on the chemical potentials, an ef- 
fective temperature dependence enters through an exter- 
nal field jl^] . Because of the flavor symmetries and anti- 
unitary symmetries of the partition function the temper- 
ature field has to satisfy the following conditions 

• The temperature term must be real. 

• It may not break the flavor symmetry. 

• Its eigenvalues are given by ±inT. In this article, 
we will consider only n = 1. 

A real temperature term satisfying these constraints is 
given by the matrix 



;(T) = 



T 

-T 



(3.2) 



For vanishing chemical potential, such a temperature 
term leads to a second order phase transition 44] . We 
wish to stress that the temperature term here is different 
from the one for QCD with three colors. Naively taking 
the same temperature dependence as for the three color 
case breaks the flavor symmetry and leads to a very dif- 
ferent phase diagram |5ll ]. 

The matrix representation for the Dirac operator in 
the random matrix model with two flavors is 
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D = 



mi 






Ml 









+ uj{T) + /ii \ 

W + i^j{T) + H2 

TOi 

"12/ 



(3.3) 



where is a real n x n matrix. The probabiUty distri- 
bution of its matrix elements is given by 



P{W) 



exp 



-G'^TtW^W 



(3.4) 



When we write the determinant as an integral over 
fermionic variables, the random matrix partition func- 
tion becomes 

Zbmt = I VWY[dij^d'4;^P{W)exp[-{pDil;]. 
■' / 

(3.5) 

We will also include source terms for a pion condensate, 



and for a diquark condensate, 



(3.6) 



(3.7) 



where the anti-symmetric Pauli matrix T2 acts in flavor 
space. 

For a fixed number of flavors {Nf = 2), the partition 
function is a function of /xi, /Lt2, T and mi, m2. A, j. While 
the masses act as source terms for the chiral condensates, 
they are also parameters of the model. 



IV. EFFECTIVE PARTITION FUNCTION 

Due to the unitary invariancc of the random matrix 
partition function, it is possible to rewrite it in terms 
of effective degrees of freedom. We will do this below, 
and then analyze the resulting effective partition func- 
tion, first to lowest order in chiral perturbation theory 
and then beyond this domain. 

Starting from the random matrix partition function 
(|3.5() . we perform the Gaussian integrations over the 
elements of the random matrix. To decouple the re- 
sulting four-fermion terms, we perform a Hubbard- 
Stratonovich-transformation and introduce an integra- 
tion over the mesonic low-energy degrees of freedom. The 
integrations over the fermionic variables can then be per- 
formed exactly. The resulting partition function in terms 
of the new variables is 



where 



Z''^ = j 2?Aexp[-£(A,At)], 



(4.1) 



C = ^G^Tr^^t _ ^logdetQ'. (4.2) 



The dependence on the chemical potentials and the tem- 
perature is contained in the 8A^/ x 8A^/ matrix Q' 



Q' 



( A'^+M ^iih+^isB -T \ 

At + M T mh+HBB 

-{fills + PLbB) -T A-fAft 

V T -{fiih + fiBB) A + Mt ) 



(4.3) 



The matrix ^ is a complex, antisymmetric 2Nf x 27V/ 
matrix. All the source terms are contained in the 27V/ x 
27V/ matrix 7lf, which for 7V/ = 2 is 



M = 



/ m —ij A 

— m A ~ij 

ij —A m 

V —A ij — m 



(4.4) 



where we have set mi = m2 = m. From here on we 
will only consider quarks with equal mass. The chemi- 



cal potentials enter as the coefficients of the matrices for 
baryon and isospin charges, 

B = diag(l, -1, 1, -1), h = diag(l, -1, -1, 1). (4.5) 

This particular form of the charge matrices is due to the 
appearance of quarks and conjugate anti-quarks in the 
spinors transforming under the enlarged SU {2Nf) sym- 
metry, which have opposite charges with respect to the 
chemical potentials. 

This effective Lagrangian can be simplified by a uni- 
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tary transformation of the off-diagonal blocks, which 
leaves the value of the determinant unchanged. Exploit- 
ing the block structure, the effective Lagrangian 14.2|l can 
be expressed as 

£ = ^G^TrA^t „ ^iogdetQ+g_, (4.6) 



where Q± are the ANf x 4Nf matrices 

A'<+M ^lBB + ^lIh±iT 



Q±- 



(4.7) 



The effective partition function with Lagrangian (|4.6|l is 
invariant under 




A 



Ji3 



-A, (4.. 



(4.9) 



This transformation thus corresponds to an interchange 
of fj.B and fii and of the diquark condensate and the pion 
condensate. 

We stress that this effective Lagrangian is an exact 
mapping of the original random matrix partition func- 
tion. 



A. Observables 

We wish to study the partition hmction and obtain the 
phase diagram with regards to four order parameters: the 
chiral condensates (uu) and {dd) for the two quark fla- 
vors, the pion condensate ^ ((M75(i) — ((J75U)) and the di- 
quark condensate ji {(d'^u*) — (u'^d*) + {(Fu) — {u^ d)). 
We identify these observables by taking the derivatives of 
the partition function with respect to the sources. After 
partial integration the expressions for the chiral conden- 
sates for the two flavors are given by 

1 



{dd) 



^a„,logZ- 

^a™, logz^*^ 

GV2. 



A 



43 



A* 



(4.10) 



For the pion condensate we find in terms of expectation 
values for the elements of the matrix A 



\{{ui-,d) 



{dlbu)) 



i^a.logZ- 



G\{A 



32 



+^32 ^23 

^ G'p 
The diquark condensate is 
1 



A=0 

A2-. + A41 - Au 

AXi — A\^ 



(4.11) 



i((d^u*) - {uU*) -f {d^u) - (u'^d)) 



1 

An 



djlogZ' 



cS 



3=0 



G^-i{Ai3 - A31 + A24, - A42 
I24 + ^42) 



^13 + ^31 



G'A. 



(4.12) 



In our convention, expectation values of the diagonal 
blocks correspond to the chiral condensates, and expec- 
tation values of the off-diagonal blocks correspond to the 
diquark and pion condensates. 



B. Chiral Lagrangian 

The chiral Lagrangian is completely determined by the 
flavor symmetries of the microscopic theory and the pat- 
tern of the chiral symmetry breaking. Thus, it must 
be in agreement with the Lagrangian of chiral pertur- 
bation theory. For vanishing chemical potential and zero 
temperature, it has been shown that the random matrix 
partition function and the static part of the chiral La- 
grangian in QCD are indeed equivalent (40,. 41. ,42]. 

To show this equivalence in the present case, we will ex- 
pand the Lagrangian H4.6() in terms of the quark mass and 
the chemical potentials around the saddle point solution 
at m = fij = ul r = 0. In accordance with the usual power 
counting rules [13 we use m = O(e^) and /u/^s = 0{e) 
and expand to second order in e. In contrast to the case 
Nc — 3, the baryon chemical potential does not drop out 
of the effective Lagrangian. The baryon chemical poten- 
tial couples to low energy degrees of freedom in form of 
diquark states, which carry a baryon charge. 

For vanishing mass and chemical potentials, the saddle 
point equation corresponding to l|4.6|l is given by 



^A = (AtA + T2)A, 
which has the nontrivial solution 



(4.13) 



A = lfT(r)i]- i(i-G2r2)i/2£^ (4.14) 

where E is any complex, antisymmetric, unitary ma- 
trix. Expanding the effective Lagrangian (|4.t)|) around 
this saddle point to order ©(e^), we find 

Tl 

'C(S) = - {ca{fIB,^J■I,T) 
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-GCT(T)Tr(SM + E^A/t) 

+0{e^)) . (4.15) 

We note that the temperature does not break any flavor 
symmetries and has no effect on the matrix structure of 
E, as we expect for a correctly implemented tempera- 
ture dependence. Its only effect is to change the overall 
magnitude of the condensates. This result agrees with 
the static part of the chiral Lagrangian derived in [l3| . 
which is (in our notation): 

£stat(S) = -:^^Tr(EM + EtMt) 



Tr(I]t(^,/3 + ^BB)E(/i//3 + MsS)). 



(4.16) 

This Lagrangian has been analyzed in detail in fl^ . In 
the remainder of this subsection we give a brief review 
of the main results of this paper. Let us first discuss 
the symmetries of the Lagrangian. In the vacuum, the 
SU{2Nf) symmetry is spontaneously broken according 
to SU{2Nf) Sp{2Nf), in the same way as it is broken 
by the mass term. At nonzero baryon chemical potential, 
the symmetry is SU{Nf) x SU{Nf) and is broken by the 
diquark condensate to Sp{Nf) x Sp{Nf). If in addition 
the mass is nonzero, the symmetry at /it ^ is C/(l) x 
SU{Nf) which is spontaneously broken by the diquark 
condensate to Sp{Nf) [lo|. 

The mass term and the terms for baryon and isospin 
chemical potential are individually minimized by E- 
matrices corresponding to a chiral condensate, a diquark 
condensate and a pion condensate, respectively. To min- 
imize the Lagrangian, we use an ansatz consisting of a 
linear combination of these individual minima, 

E — cosaEo- + sina(cos77EA + sinr^Ep), (4.17) 

which is again antisymmetric and unitary. These matri- 
ces, which correspond to the source terms introduced in 
(|^ . l|X7|) . and are given by 
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(4.18) 



Although the baryon chemical potential term is invariant 
under the J7b(1) rotations generated by B, this symme- 
try is spontaneously broken by the diquark condensate. 



violating baryon number conservation and giving rise to a 
massless Goldstone boson . In the same way, the pion 
condensate breaks the I3 generated symmetry, which also 
leads to a violation of isospin number conservation and 
a Goldstone boson (3.15^. 

Using the ansatz for E, the effective Lagrangian 14. 6|) 
becomes 

C = nNf [1 + logG^ - 2mGcosa 

-|-G^(cos^ ai^j + n%) - sin^ acos(2?7)(^| - p.])) 
+ 0(e2)] . (4.19) 

When this is minimized with respect to 77, for a ^ Q and 
11% ^ ^j, we find 



T] ~ for ^5 > fij 
V= J for ^1 < fij, 



(4.20) 



where the first case corresponds to the presence of a di- 
quark condensate, and the second to that of a pion con- 
densate. Thus, there is a first order transition between 
these two phases on the lines = Minimizing now 
with respect to a, we find 



1 



2G fi- 

m 1 
2G^ 



for 

,,2 



9 m , , m 
• <2g""^^^<2G 



2G 

1 r 2 ^ 1 2 2 

^ for ^J■B> 2Q a-iid -"s > 

^ and fi% < 11]. (4.21) 



for /zj > 



2G 



Since the critical chemical potential for the condensation 
of Goldstone bosons in our convention is half the pion 
mass, /ic — m7r/2, we can identify to this order in the 
expansion parameter 



(4.22) 



From the analysis of the chiral Lagrangian we thus 
have obtained the following picture of the phase diagram 
at zero temperature: below a chemical potential of half 
the pion mass, the chiral condensate is constant. Above 
this critical chemical potential, depending on the relative 
magnitude of isospin and baryon chemical potential, pion 
or diquark condensation sets in and the chiral condensate 
drops quickly. The phases of pion and diquark conden- 
sation phases are separated by first order transition lines 
at Aif = 



C. Effective Potential 

We now wish to analyze the random matrix model to 
all orders in the mass and chemical potentials. To this 
end, we make an ansatz for the matrix A, and obtain an 
effective potential for the condensates in our model. As 
in the case for three colors, the ansatz for A is guided 
by two observations: From chiral perturbation theory. 
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we know that above a critical chemical potential pion 
and diquark states will condense. We also expect that 
the presence of two flavor dependent chemical potentials, 
which explicitly breaks the flavor symmetry, will lead to 
a different behavior of the chiral condensates for the two 
flavors. This leads to the following ansatz for A, 






-CTl 


-iA 


-P 


c^i 





-P 


-iA 


iA 


P 







P 


iA 


0-2 






(4.23) 



With this ansatz, evaluating the determinant we obtain 
the effective potential for the condensates: 



-C = G^{crj+al + 2A^ + 2p^)~l-TT\ogQ+Q^, (4.24) 
n 2 

where 

detQi = {[{ai+m + Hi±iT){a2+m~ iJ,2TiT) + + A'^] 

X [((Ti +m- fii± iT){a2 + m + iT) + + A^] + AA^ 1^11x2} 

{[((Ti + m - /ii =F iT){a2 + rn + ^2 ± iT) + + A^] 

X [((Ti + m + ^1 =F iT){(72 + m- p2±iT) + + A^] + 4AV1M2} , (4.25) 



Although det Q+ and det Q- are complex conjugate to 
each other by construction, they turn out to be real them- 
selves. With the ansatz (|4.23|l the symmetry (|4.8|) corre- 
sponds to the interchange (for j = A = 0) 



A 



-A, 



P2 



-P2- 



(4.26) 



The effective potential (|4.24ll is indeed invariant under 
this symmetry, as well as the interchange of the flavor 
indices, and fif ^ —pf for both f — 1,2. These are the 
symmetries of the QCD partition function, as shown in 
section ^1 They allow us to restrict our study to one 
half quadrant of the /XB-/i/-plane. The other sectors of 
the phase diagram are easily obtained by applying these 
symmetries. The transition between the phases with di- 
quark and pion condensation must necessarily be of first 
order and take place on the lines /i| = p^ . 



V. PHASE DIAGRAM 

In this section, we will determine the complete phase 
diagram of the random matrix model in the space of the 
chemical potentials, ps and pi, and the temperature T. 
To this end, we solve the saddle point equations 

t-' f-' M-o <-) 

and determine the expectation values of all the conden- 
sates in the saddle point approximation. For vanishing 
temperature, we will obtain analytical solutions in the 
chiral limit as well as for finite quark mass. For nonzero 



temperature, we will present analytical results only in 
the chiral limit and study the potential numerically for 
finite quark mass. 



A. Chiral limit at zero temperature 

First, we investigate the phase diagram in the pi-p2 
chemical potential plane at temperature T = 0. In 
the chiral limit, the Goldstone bosons arising from the 
spontaneous breaking of the chiral symmetry are exactly 
massless. Since all Goldstone bosons carry a charge that 
couples to one of the chemical potentials, this has a di- 
rect bearing on the phase diagram. At zero temperature, 
the chiral condensate will be rotated completely into a 
bosonic condensate as soon as one of the chemical poten- 
tials becomes nonzero. 

a. Chiral condensates and chiral restoration transi- 
tion. In the phase where the bosonic condensates van- 
ish, A = p = 0, the effective potential greatly simplifies 
and becomes a sum over the contributions of the chiral 
condensates for the two flavors. 



I 



-C 



E 

/=1,2 



G 



l0g(CTj 



2 \2 



(5.2) 



Exactly as in the case of three color QCD, the saddle 
point equations 

af[G^ia^f-p})-l] = 0, / = I,2, (5.3) 

have the solutions 

af = 0, / = 1,2, (5.4) 



8 



'^Z = G^+^/' / = 1'2, (5.5) 

which correspond to a phase with restored and broken 
chiral symmetry, in this order. The contributions to the 
free energy for one flavor are 

Qfifif) = -\og^i} (5.6) 
17/(/i/) = l + logG'+fi}G^ (5.7) 

for each of these phases, respectively. The chiral symme- 
try restoration phase transitions take place at the chem- 
ical potential for which the two free energies coincide. It 
is given by the solution of 

1 + fi}G^ + log fijG^ = 0. (5.8) 

We note that this condition depends on only one of the 
two chemical potentials and is strictly constant in the 
other. The behavior of the two flavors in these phases is 
completely independent. This changes only in the pres- 
ence of either a pion or diquark condensate, which couples 
both flavors. As a result, the phases with one nonzero chi- 
ral condensate extend along the chemical potential axes. 
Only in the center of the phase diagram, where these 
strips overlap, could a phase exist for which both chiral 
condensates arc nonzero. 

b. Phases with bosonic condensates. At the point 
fj-i — fJ-2 — 0, the solution of the saddle point equation 
is degenerate. It satisfies (TiO'2 + + — l/G^, with- 
out a preference for a specific direction in the space of 
the condensates. Due to the absence of flavor symmetry 
breaking terms we expect txi = cr2 • This is the only point 
where all condensates can be simultaneously nonzero. 

In the chiral limit, the chiral condensates vanish at 
nonzero chemical potential in the presence of a bosonic 
condensate. As we argued above, the phases with di- 
quark and pion condensation must be separated by a flrst 
order transition. This is confirmed by the saddle point 
equation: the only solution for which A 7^ and p 7^ 
requires that either — or fi2 — 0, which is equivalent 
to /ij — ^\ and thus coincides with the expected first 
order transition line. Therefore, we always have either a 
pion or a diquark condensate, and the effective potential 
in these two cases is reduced to 

-£ = G^2p' -Xogip" - iiui2)\ (5.9) 
n 

-C ^ - log(A2 ^1/12)^ (5.10) 

n 

The corresponding saddle point equations become 

GW-M1M2)-P = 0, (5.11) 
G2A(A2 + ^i/i2)- A = 0, (5.12) 

and have the respective solutions 

P =0, p2 =_L+^^^2 (5.13) 
A =0, A2 =i^-;,i^2. (5.14) 




jjLlG 



FIG. 1: Phase diagram in the ^1-^2-plane at T = and van- 
ishing quark mass mG — 0. First (second) order transitions 
are indicated by solid (dashed) lines. Phases are labeled by 
the nonvanishing condensates. All condensates vanish in the 
corner regions, where chiral symmetry is restored. 

Second order phase transition lines are given by the con- 
ditions /X1/Z2 = — 1/G^ and = 1/G^. The free ener- 
gies for the two solutions are 

rip = 2(l-HlogG2+/ii^2G2), (5.15) 
= 2{l + \ogG^ ~ fiifi2G^). (5.16) 

By comparing the free energies, we observe that in the 
quadrants with fj,ifj,2 > 0, the diquark condensate is 
nonzero, and that for /ii/L(2 < the solution with a pion 
condensate has the lower free energy. Equating the free 
energies, we confirm again that first order transition lines 
are given by fiifi2 = 0. 

c. Phase diagram. To find the complete phase di- 
agram, we use the results we obtained directly for the 
second order phase transitions, and we compare the free 
energies of all phases to determine the first order tran- 
sitions. We conclude that, for massless quarks, phases 
with condensed bosons are always favored above a phase 
where both chiral condensates are non-vanishing. In ad- 
dition to the first order transitions between the pion and 
diquark condensate phases, we find first order transitions 
from these to the phases where only one of the chiral con- 
densates is different from zero. 

The phase diagram is shown in Fig. ^ At the origin, 
the solution to the saddle point equation is degenerate, 
and only the magnitude of the condensates is fixed. From 
the origin, a pion condensation phase extends into the 
quadrants where /ii/i2 < (/Uj > and a diquark 

condensation phase extends into the quadrants where 
> {fij < ^i^g). The domain of these phases is 
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limited by a line of second order transitions (see dashed 
curve in Fig. which in each quadrant is given by the 
solutions of 



1±GViA*2 =0. 



(5.17) 



These transitions can be interpreted as saturation tran- 
sitions which also occur in lattice QCD because of the 
finite number of available fermionic states . Our sat- 
uration transition is somewhat different. Because we are 
working in the thermodynamic limit, the number of avail- 
able states is infinite as well. However, since the level 
spacing of the relevant matrices ~ 1/n, the support of 
the spectrum remains compact. It is the non-analyticity 
that occurs because of the extreme edges of the spectrum 
that gives rise to the saturation transition. The phases 
with broken chiral symmetry and only one of the chiral 
condensate nonzero appear along both chemical potential 
axes. They are separated from the Goldstone condensa- 
tion phases by first order transition lines. 

The phase diagram exhibits a fourfold symmetry, along 
both the fiB- and the /x/- axis. This is in contrast to the 
twofold symmetry with respect to the ^s-axis that we 
found in the phase diagram of three color QCD, where 
only a pion condensate appears [s^ . Fundamentally, this 
is due to the flavor symmetry in the two color case which 
is enlarged beyond the symmetry in the three color case. 
However, the upper- left quadrant in the /is-M/~plaiie is 
strictly identical in our random matrix model for both 
two and three colors |33| . The important point is that 
both chemical potentials couple to Goldstone bosons of 
the theory. Therefore, the two chemical potential terms 
in the present case are on equal footing. They lead to 
two different possible symmetry breaking patterns. For 
three colors, there is only one chemical potential term, 
since only the isospin chemical potential couples to the 
Goldstone modes, and thus there is only one way to break 
the flavor symmetry. 



B. Nonzero quark mass at T = 

The phase diagram in the presence of a flnite quark dif- 
fers from the one obtained in the chiral limit. As in the 
case of three colors, this is due to two effects: First, the 
chiral condensate becomes nonzero because of the explicit 
breaking of chiral symmetry. Because of this, phase tran- 
sitions may turn into crossover transition. Second, previ- 
ously exactly massless Goldstone bosons become massive. 
This primarily affects the low energy regions of the phase 
diagram. In addition, the Goldstone condensation phases 
become qualitatively different and are now mixed phases 
with nonzero chiral condensates. In spite of this, at zero 
temperature we can still solve the saddle point equations 
exactly, and we will give the results below. 

a. Chiral condensation phases. For a small quark 
mass, the behavior of the chiral condensate is still domi- 
nated by the spontaneous breaking of the chiral symme- 
try. Thus, a natural way to treat the chiral condensates 



in this case is to expand in the quark mass about the 
solutions for m = 0. 

The effective potential that describes the phases with- 
out bosonic condensates is given by 



^ J2 GV^-ilog[(a^ +m) 



2 2l2 



/=1,2 



(5.18) 



For each flavor we have to solve the saddle point equation 

G'^afiiaf + jnf - fi}] - {<7f + m) = 0. (5.19) 

As expected, the solution af — no longer exists, and 
we find by expanding in m (without assumptions about 
the magnitude of fif) the new solution 



G2 



TO + 0{m^). 



(5.20) 



Expanding in mG in the phase with spontaneously bro- 
ken chiral symmetry around the second solution with 
G^cTf — 0{G), we find the solutions 



Gaf = ±.^l + njC 
+0{m^G^) 



1 



2{l + li}G^) 



- 1 mG 



(5.21) 



The respective free energies for one flavor of these two 
phases with spontaneously broken and restored chiral 
symmetry are to lowest order in to given by 



% = 1 + tijG' +logG^T2mG^l + tijG^l5.22) 
flf ^ -lognj. (5.23) 

If these two phases are not separated by a Bose-condensa- 
tion phase and meet directly, the chiral restoration tran- 
sition is of first order. The chemical potential at the 
transition is obtained by matching the free energies of 
the two phases. We find 

mG 



Ai/,cG = fJ-cG 1 



1 + M?G2 



+ 0{m^G') ,(5.24) 



where /XcG is the critical chemical potential in the chi- 
ral limit determined from H5.8|l . The shift in the chiral 
restoration transition is proportional to the quark mass 

TO. 

b. Goldstone boson condensation phases. For finite 
quark mass, the critical chemical potential for condensa- 
tion of Goldstone bosons becomes nonzero due to their 
finite mass. The symmetry with regards to the chemical 
potentials fj,i and fis remains unaffected. In particular, 
above the critical chemical potential, we still expect to 
see either a phase with diquark condensation or one with 
pion condensation depending on the relative magnitude 
of these two potentials, and a first order transition sep- 
arating them. Since the results for both Goldstone con- 
densation phases are related by this symmetry and are 
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thus very similar, we will give explicit results for one half 
quadrant only. 

In both phases, due to the finite quark mass, the chiral 
condensates do not vanish. Therefore, we have to solve 
three saddle point equations simultaneously, for the chi- 
ral condensates of both flavors and for either p or A. In 
the region in which > {fisl, we thus have A = 0. 
The solutions of the saddle point equations are given by: 



m 1 



1 



0-1 



cr2 



2 GV? - 



1 . , sMs 



m 1 



1 



fJ-B 



1 



2 2,2 2Ms 

G M7 



(5.25) 



From the relation between the chiral and the pion con- 
densate, we find the explicit result for p, 



2 2 2,2 

P ^ Pb- ^-1+ Pi 



1 



1 



m"' 1 



1 



2 Pb 1^/2 2 , 



2 4 G4 {fij - m2)2 

11 1 



(5.26) 



The difference between the two chiral condensates de- 
pends on the ratio of /zg and fj,j: 



^ Pb 
PI ' 



The free energy of this phase is given by 



np{m,HB,P-i) = 2(l + logG2 



-m^G^^ 
Pi 



+ {^Jil - m?)G2 
log 



P'l 



P] 



(5.27) 



^G^) 



(5.28) 



The condition — q describes a line of second order 
phase transitions at which the bosonic condensate van- 
ishes. For finite quark mass, this condition has more than 
one solution. For fiB — ^, the two solutions for p?j are 

p]cG' - '^ + lm'G^+0{m'G') (5.29) 
^l],G^ = l + -^m'G^ + 0{m^G^). (5.30) 

We note that the solution for which fij^G^ ~ mG/2 + 
0{m^G^) agrees to leading order with the result which 
was obtained from chiral perturbation theory JX}\ . For 
chemical potentials below the critical value of at 
which diquark condensation occurs, the pion mass is ex- 
pected to be independent from /xs 13] . Thus, the critical 
value for fij should be exactly constant as a function of 
IJ,B- In our random matrix model we find that fij c is 
only approximately constant as a function of /is. This 
is due to the fact that, in random matrix models, the 
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FIG. 2: Phase diagram in the /ii-/i2-plane at T = for a 
finite quark mass mG = 0.1. First (second) order transitions 
are indicated by solid (dashed) lines. The difi'erent phases are 
labeled by the non-vanishing condensates or the condensates 
which have large expectation values. The finite quark mass 
leads to a phase with chiral condensation only in the center 
of the phase diagram. In the corner regions, both the diquark 
condensate and the pion condensate vanish, and the chiral 
condensates are of 0{m). 



chemical potential affects the condensates already below 
the threshold, even at zero temperature |R l46l Isif. In 



15 1| it was thus argued that the random matrix result 
should be interpreted as the critical contribution to the 
free energy, on top of a noncritical part not contained in 
the model. 

For small baryon chemical potential, small isospin 
chemical potential and quark mass, the chiral conden- 
sates in the pion condensation phase become 



(Tl = (72 



2 I 2 
(T + p 



1 

CP 



+ 0{m) (5.31) 
0{m), (5.32) 



2G2p2 



which agrees with the results H4.21|l obtained in chiral 
perturbation theory. In this region of the phase diagram 
the effects of a finite quark mass are most significant. For 
large values of the chemical potentials both the pion and 
diquark condensates vanish for any value of the quark 
mass. The second order transition to this phase is to 
leading order independent of the quark mass. 

c. Phase diagram. As in the chiral limit, in order 
to obtain the complete phase diagram we combine the 
results for the second order transitions with an analysis 
of the free energies of the different phases. The result- 
ing phase diagram in the /ii-/i2-plane is shown in Fig. [21 
for a quark mass of mG = 0.1. The four- fold symmetry 
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FIG. 3: Phase diagram in the /ii-/i2-plane at T = for a 
finite quark mass mG = 0.2. First (second) order transitions 
are indicated by solid (dashed) lines. The different phases 
are labeled by the condensates which are non-vanishing or 
have a large expectation value. As described in the text, for 
a quark mass larger than mG ~ 0.184 the pion and diquark 
condensation phases no longer have a boundary in common. 
In contrast to the situation shown in Fig.|21 the chiral first or- 
der transition lines therefore appear also if one of the chemical 
potentials /i/ is small. 

remains the same as in the case of mG ~ 0. In particu- 
lar, the transitions between the phases with the different 
Goldstone boson condensates are still of first order. The 
main difference from the chiral limit is the appearance of 
a phase in the center of the phase diagram in which the 
chiral condensates for both flavors are nonzero. This is 
a consequence of the finite critical chemical potential for 
the condensation of bosons with a finite mass. 

A second effect of the finite quark masses is that the 
chiral restoration transitions outside the regions of Bose 
condensation are shifted slightly from their positions in 
the chiral limit. As can be seen from H5.24|l . the correc- 
tion is of order mG. Thus, while changing the nature 
of the phase transition considerably in the regions of low 
chemical potential in the center of the phase diagram, a 
finite quark mass has only a small effect on the phase 
diagram at large chemical potential. 

For leading order in the mass expansion, we find that 
the first order transition between the phases with p ^ 
and A ^ disappears at the point 1 + log2mG — 
0. Therefore, for masses larger than mG ~ 0.184, the 
pion and diquark condensation phases are not contiguous 
anymore. In this case, if either fii or ^2 is kept small 
and the chemical potential for the other quark flavor is 
increased, one encounters the chiral restoration transition 
without passing through cither the pion or the diquark 



condensation phase. 

C. Nonzero T at m = 

1. The Hi = 0- and fiB = 0-planes 

Since the presence of a temperature term does not 
change the symmetries, even at finite temperature the 
two phases with the different Goldstone condensates are 
still related by an interchange of the two chemical po- 
tentials. We can again treat both the pion and diquark 
condensation phases in parallel. In either the fii — 
or the = plane, the presence of only one of these 
chemical potentials restores the symmetry between the 
two quark flavors. Consequently, there is only one chiral 
condensate a = ai = 1T2 ■ The results for the chiral phase 
transition were first obtained for /i/ = in The re- 
sults in the plane = coincide with those for Nc = 3 
at nonzero isospin chemical potential 33] where the same 
effective potential was found. Therefore, and mainly for 
completeness, we briefly review the main results. The 
result for the diquark condensate can be obtained by the 
simple substitution /i/ //_b and p ^ — A. 

In the chiral limit the critical chemical potential for 
pion condensation of this model is independent of the 
temperature. For any nonzero value \pi\ > 0, the system 
will be immediately in the pion condensation phase (via a 
first order transition from the degenerate point pi = 0). 
The effective potential for this phase is 

-£ - - log[(p2 +pj+ TY] , (5.33) 

n 

and the saddle point equations are solved by 

- ^~f^'i-T^ (5.34) 

p = 0. (5.35) 

The phase diagram is now simple: A second order tran- 
sition fine is given by the circle G'^{pj + T^) = 1. The 
free energy of this phase with cti = (T2 — 0, given by 

Qp = 2[1 + logG^ ~G^{pj + T^)], (5.36) 

is always lower than that of the chiral broken phase, and 
thus in the chiral limit, the chiral condensates do not ap- 
pear in this plane (see the upper panel of Fig.0J. There- 
fore, the tricritical point of the chiral restoration transi- 
tion found in the ^s-T-plane at pi = for iVc = 3 Q 
does not appear at all in the /j,/-T-plane given hy pB = 0. 
Using the symmetry of the phase diagram, we see that, 
in contrast to the case Nc = 3, for Nc = 2 this tricriti- 
cal point is not present in the /is-T-plane with pj = 0, 
either. 

2. The complete pi -pB -plane at T ^ 

In the chiral limit, we can still solve the saddle point 
equations analytically. Since the first order phase transi- 
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tion between the pion and diquark condensation phases 
is not affected by the temperature, we continue to have 
phases of pure diquark and pure pion condensation which 
meet at the first order transition hues fii — or fi2 — 
(i.e. /i/ = zb/is). There are no phases in which a pion 
condensate or a diquark condensate coexists with a chiral 
condensate of any flavor. 

a. Chiral condensation phases. The phase with only 
chiral condensation occurs along strips parallel to the 
^1 = or ^2 = axis (see Fig.0. In these phases only 
one of the chiral condensates is non- vanishing. The effec- 
tive potential simplifies to a sum over two independent 
contributions of both flavors, 



rp2 



(5.37) 



and coincides with the effective potential found for QCD 
with three colors. For each of the two flavors, the chi- 
ral condensate depends only on the temperature and the 
chemical potential for that flavor. In the chiral limit, 
the magnitude of the chiral condensate is given by the 
solution of the saddle point equation 



The tricritical point, where the phase transition line 
changes from first to second order, appears where the 
coefficients of u/ and aj in H5.38|l vanish. We find 



4 ' 
V2 + 1 



(5.44) 



Since these results are independent of the chemical poten- 
tial for the other flavor, all transition lines are straight 
lines in the /ii-/i2-plane at fixed temperature. In fact, 
the transition lines for one flavor in this plane have to 
be perpendicular to the corresponding lines of the other 
flavor. 

b. Goldstone boson condensation phases. Once 
again relying on the symmetry of the effective potential 
(|4.24|) . it is sufficient to analyze the pion condensate. 
Results for the diquark condensate simply follow from 
interchanging fii and fis- The effective potential in the 
pion condensation phase at finite temperature is given 
by 



n 



2GV^ 
log[(p2- 



ry + AfijjT^]. (5.45) 



G2 



Its solutions are given by 
af - 



0. (5.38) 



(5.39) 



1 



'f 



2G2 



(5.40) 



and due to the solution af = a first order transition to 
the phase with restored chiral symmetry is possible. For 
any one flavor / = 1,2, this flrst order transition takes 
place where the free energies of the two solutions, given 
by 



% = -logV^+T^) 



(5.41) 



log G' + G' {^ij -T')±-Jl- {AG^fifT) 



-,^o,[l±yi-i,G^,fTr 



(5.42) 



respectively, become equal. The solution with 

the negative square root is not a global minimum of the 
free energy, and is in fact never realized. We find that 
the solution ct/ ^ approaches the solution tr/ = at 
the second order phase transition line in the /^/-T-plane 
given by 



(5.43) 



To explore the structure of the boundaries of the domain 
of this phase, we analyze the saddle point equation 



K2G2 



2 



4M|r2 



= 0. 

(5.46) 



Besides the trivial solution p = 0, the saddle point equa- 
tion has the solutions 



= ^ + - M? - ± ^ Vi - (4G2^Br)2 



(5.47) 

of which only the one with the positive sign in front of the 
square root is a minimum of the free energy. A second 
order transition line is given by the condition that this 
solution coincides with the trivial solution, 



+^(Ml-/i?-r2) = o. 



(5.48) 



There is always the possibility of a first order transition 
to a phase with p = 0. Indeed, for a certain range of 
parameters, where the coefficients of both p^ and p in eq. 
(|5.46|) vanish, there exists a tricritical point at which the 
transition changes order. But when the complete phase 
diagram is considered, including all other condensates, 
the solutions that lead to this tricritical point do not ap- 
pear: at this point the free energy of one of the chiral 
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condensation phases is lower than that of the pion con- 
densation phase, and a first order transition takes place 
before the point is reached. The free energy in the pion 
condensation phase is given by 



1 



2^-+logG^ + G^(/4-A^f-T^ 



■log 



(5.49) 



and we determine the complete phase diagram by com- 
paring the free energies of the different phases. Using the 
symmetry (|4.8|l . all the above results also apply to the 
diquark condensate. 

c. Phase diagram. At finite temperature, the phase 
diagram in the /Ui-/i2 plane still has the fourfold symme- 
try as seen at T = in Fig. ^ As can be seen from the 
second order transition condition H5.48|l . the boundaries 
of the regions with diquark condensation or pion conden- 
sation move gradually towards the origin fis = fJ-i = 
with increasing temperature, and shrink rapidly to a 
point when TG approaches TG = 1. The regions with 
a nonzero chiral condensate also shrink gradually with 
increasing temperature. When the temperature reaches 
the tricritical point given by eq. (|5.44() . the transition 
to the chirally restored region becomes second order. As 
the temperature approaches TG — 1 this region rapidly 
collapses onto the axes /xi = and ^2 = 0. All conden- 
sates vanish and the symmetries are completely restored 
for TG > 1. 



D. Finite quark mass at nonzero T 
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In general, for m 7^ the saddle point equations cannot 
be solved analytically when iij, fis said T are all different 
from zero. However, as outlined in the previous section, 
the symmetries of the partition function under an inter- 
change of the chemical potentials for the two flavors, or a 
change in their sign, are unaffected by the presence of the 
temperature term. Thus, the four symmetry axes /x/ = 0, 
fiB = and /i/ = ±/iB remain at finite temperature, and 
the first order transitions between the condensates A and 
p at fii = ±fiB remain. The structure of the phase di- 
agram in this case can be deduced from the analytical 
results at T = and from the planes where one of the 
chemical potentials vanishes. 

1. The fiB ~ 0- plane 

Although an analytical solution of the saddle point 
equations is not possible for both chemical potentials and 
the temperature nonzero, in order to have a complete pic- 
ture of the phase diagram it is instructive to consider the 



FIG. 4: Phase diagram in /i/ and T in the /j-b ~ 0-plane, both 
for the chiral limit (top) and a finite quark mass mG = 0.1 
(bottom). First (second) order transitions are indicated by 
solid (dashed) lines, the dotted lines in the case of finite quark 
mass mark the crossover. 



plane hb =0. Here, of course, the chiral condensates 
for both flavors coincide again. This region of thephase 
diagram has been studied at zero temperature in [j| and 
by means of a chiral Lagrangian for small values of 
fij. Finite temperature effects have been included in a 
one-loop calculation in chiral perturbation theory in [l^ 
and in the random gauge model of |5lj |. 

a. Chiral condensate. Below the critical value of /i/, 
chiral symmetry is broken by the chiral condensate, and 
at high temperature it is restored after a crossover tran- 
sition. For m ^ 0, it is most transparent to solve the 
saddle point equations for the chiral condensate by ex- 
panding in the quark mass m. In the broken phase, the 
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FIG. 5: Phase diagram in the /ni-/i2-plane at fixed values of TG and for a finite quaxk mass of mG = 0.1. Values for the 
temperature TG are given in the panels. First (second) order transitions are indicated by solid (dashed) lines, the dotted lines 

indicate a crossover. Phases are labeled by the non-vanishing condensates, rcsp. the condensates with large expectation values. 
In the corner regions, the pion and diquark condensates vanish, and both chiral condensates become small. 
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chiral condensate a = ai = a2 is given by 

1 a^-fij+T^ -SG^T^nj 



a — ao + m 



2al 1 - (4G>/r)2 



+ ^\/l-(4GW)^. (5.51) 

m + 0{'m ) increasing temperature, the chiral condensate de- 

creases. In the high temperature phase, the chiral con- 
densate can be expanded about the trivial solution (5.39) 



where cro is the solution of the saddle point equation for as 
m = 0, 

^2 _ 1 I ,,2 T-.2 

•^0 - 



a = —m 



G^{lij+T^)^ + Hj-T^ 



+ 0{m^). (5.52) 



16 



b. Pion condensate. As in the case T — 0, for fi- 
nite quark mass the chiral condensate does not vanish in 
the pion condensation phase. In this mixed phase, we 
can solve the saddle point equations analytically. The 
solution given by 



2G2 



2 I 2 

a + p 



1 

-Q2 



m 



(5.53) 
(5.54) 



was first obtained for diquark condensation [5l|. We no- 
tice that the chiral condensate is completely independent 
of the temperature, and only the total magnitude + p^ 
of the condensates decreases with increasing tempera- 
ture. At zero temperature these equations describe the 
rotation of a chiral condensate into a pion condensate for 
increasing /i/. The explicit solution for the pion conden- 
sate obtained from H5.54|l is given by 

2 _ 1 



p?j — -m? 



4G4 {lij - m 



2\2 



Ai?-r2. (5.55) 



The phase with p ^ is bounded by a second order line 
at which p vanishes: 

2,' 2 2\ -'^ 2 

Pi{Hi -m )- -m 

-{li]+T'')[n]~m^fG^ = 0. (5.56) 

An important observation is that in this for 
Nc = 3, the boundary of the pion condensation phase 
is always a second order transition, and there is no tri- 
critical point. The complete phase diagram for p,B — ^ 
is pictured in Fig. 0] In the /i/-T-plane the critical end- 
point of the first order chiral restoration transition does 
not appear, because its position is in the region where 
the pion condensation phase has a smaller free ene rgy 
than the chiral condensation phase. We have shown |33 
that in this plane even for large quark masses, the crit- 
ical endpoint remains masked by the pion condensation 
phase. 



2. The complete hi-^lb -plane at T ^ 

To obtain the full picture of the phase diagram at 
nonzero values for the two chemical potentials, the tem- 
perature and the quark mass, we solve the saddle point 
equations numerically. The results for the phase dia- 
gram are presented in Fig. [S] It is in complete agreement 
with the expectations from our discussion of the phase 
diagram at either m = or T = 0: The finite quark 
mass mainly affects the phase diagram at low values of 
the chemical potentials, and the regions with either a di- 
quark or a pion condensate shrink with increasing tem- 
perature and finally vanish for a temperature of order 
one, TG = 0(1). 



We wish to draw the attention to one particular de- 
tail of the phase diagram in our model. From Fig. 
the diquark and pion condensation phases are no longer 
contiguous at the lines fii — and fi2 — 0, and there 
is no first order transition between them. Instead, this 
region of the phase diagram is occupied once again by 
the chirally broken phases, and the chiral phase transi- 
tion line reemerges. While the first order chiral transition 
line does not appear for /x^ = (and by symmetry not 
for fij = 0, either), we thus find that for both p,B and 
fii different from zero and an appropriate value of the 
quark mass this first order line and its critical endpoint 
can be observed. The phase diagram in the /is-T plane 
at an isospin chemical potential that is chosen to avoid 
the superfluid phases at low fis shares striking similari- 
ties with the phase diagram of QCD with three colors at 
nonzero baryon and isospin chemical potentials. 



3. The fiB-T-plane at constant pi 

To study the chiral phase transitions for both ps and 
Pi nonzero in some more detail, we turn to the phase 
diagram in planes described by constant isospin chemical 
potential. Fig. |B| shows the phase diagram in these ps- 
T-planes for fixed values of pj. The corresponding cuts 
through the phase diagram in the three parameters pi, 
Pb and T are perpendicular to the cuts at constant TG 
in Fig.O Since pi = Pb + Pi and P2 = Pb — Pi, planes 
at constant pi intersect the planes with constant T along 
lines parallel to the diagonal from the lower left-hand to 
the upper right-hand corners of the panels in Fig. |S1 

In our analysis of the plane pB = and in Fig.0]above, 
we have seen that the first order chiral phase transition 
and its critical endpoint are not present for any value of 
the quark mass, and by virtue of the symmetries this also 
holds for the plane given hy pi = 0, shown in the first 
panel of Fig. El But with increasing value of the isospin 
chemical potential the phase diagram in the pB-T-plane 
changes considerably: first, the first order chiral transi- 
tion for CTi ^ (uu) at large temperature and its critical 
endpoint reappear. With a further increase of pi, the 
chiral transition for a2 at large values of pB appears as 
well (see the second and following panels of Fig. 0) . We 
have to caution, though, that for the second chiral con- 
densate (72 this is possible only because of the saturation 
transition of the diquark condensate, which exists on the 
lattice as well as in our random matrix model. If pi is 
made larger still, the critical endpoint of the transition 
for CTi is once again swallowed by the pion condensation 
phase (see last panel of Fig.El where this is about to take 
place) . 

For a quark mass mG > 0.184 and at any temper- 
ature, the diquark and pion superfluid phases are not 
contiguous. Therefore for a sufficiently large quark mass 
(mG > 0.184) and an isospin chemical potential that is 
small enough so that the pion superfluid phase is never 
reached and large enough so that the diquark superfluid 
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FIG. 7: Phase diagram in the ^s-T-plane at fixed values of 
the isospin chemical potential fj,iG = 0.3 for a finite quark 
mass of mG = 0.2. First (second) order transitions are in- 
dicated by solid (dashed) lines, the dotted lines indicate a 
crossover. Phases are labeled by the non-vanishing conden- 
sates, or the condensates with large expectation values. This 
phase diagram is very similar for either QCD with two colors 
or for QCD with three colors. 



phase is reached only at large /is, the structure of the 
phase diagram is as in Fig. 13 The resulting phase dia- 
gram is strikingly similar to the phase diagram of the 
Random Matrix Model for QCD with three color at 
nonzero baryon and isospin chemical potentials js^ . 



VI. CONCLUSION 

In this article, we have introduced a random matrix 
model with the symmetries of QCD with two colors 
and two flavors and nonzero baryon chemical potential, 
isospin chemical potential and temperature. Contrary to 
work elsewhere in the literature, our temperature term 
respects the flavor symmetry of the partition function as 
it should. 

The random matrix partition function has been rewrit- 
ten exactly in terms of an effective Lagrangian. To lowest 
order in the chiral perturbation expansion, the effective 
Lagrangian agrees with the one derived in chiral pertur- 
bation theory 0, ^| . We have analyzed the phase 
diagram of this model for values of the mass and chemi- 
cal potential outside of the domain of chiral perturbation 
theory. It shows a rich phase structure with phases char- 
acterized by the diquark condensate, the pion condensate 
and the chiral condensates, {uu) and {dd). In the pres- 
ence of two chemical potentials, the chiral condensates 



are not necessarily equal. For a nonzero quark mass, 
in regions without pion or baryonic condensate at small 
values for the chemical potentials, the restoration of the 
spontaneously broken symmetry takes place via two sep- 
arate crossovers, across which the chiral condensates for 
the two flavors rapidly change in value. The vanishing of 
the superfluid phases at large values of the relevant chem- 
ical potentials, which is also observed on the lattice, has 
been explained as a saturation effect. On the lattice, this 
is due to the flnite number of lattice sites [23, and for 
random matrix models, this is due the compactness of 
the Dirac operator. 

The tricritical point of the diquark and pion con- 
densation transitions observed in lattice calculations 
[IE HE HE EE| is absent in the random matrix models. 
Since the Random Matrix partition function is equiva- 
lent to a mean fleld description of the phase transition, 
this is not entirely surprising. Symmetries alone do not 
determine the presence of a tricritical point. In chiral 
perturbation theory, a tricritical point is obtained from 
a one-loop effective potential (12ii. i32] . Ultimately, the 
question of the structure of the phase diagram in this 
regard must be answered by lattice calculations j^El- 

Finally, we would like to compare our results with the 
random matrix study of the phase diagram for three- 
color QCD at nonzero temperature, baryon and isospin 
chemical potentials 33]. First we notice that the upper- 
left quadrant of the phase diagram in the /is-^j-plane 
is identical in both cases. This is due to the fact that 
in both cases, the isospin chemical potential is related to 
Goldstone bosons of the microscopic theory. The upper- 
right quadrants are obviously different in the two theo- 
ries, since there is no Goldstone boson that carries baryon 
number in three-color QCD. However, the phase diagram 
for three-color QCD with three quark flavors of equal 
masses at nonzero isospin and strange chemical potential 
should be very similar to the one we have presented here. 

Second, and more importantly, we notice that it is 
possible to flnd a quark mass and an isospin chemical 
potential for which the diquark and pion condensation 
phases appear only in certain regions of the phase dia- 
gram. Namely, they are present only at low temperature 
and either very small or very large values of the baryon 
chemical potential. As shown in Figs. IHl and [7| we see 
that the flrst order transition lines of the chiral transi- 
tions and their critical endpoints are not at all affected by 
the phases with bosonic condensates. Indeed, the tran- 
sition lines are identical to those that are found in the 
random matrix model for three-color QCD. This obser- 
vation implies that the critical endpoint of three-color 
QCD can be studied in two-color QCD. For an appropri- 
ate choice of quark mass and isospin chemical potential, 
the phase diagram in the ^b-T plane are almost iden- 
tical for QCD with two colors and for QCD with three 
colors. In particular, the doubling of the critical lines due 
to a nonzero isospin chemical potential can be studied in 
QCD with two colors. 
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